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ABSTRACT 

We present the results of a numerical magnetohydrodynamic simulation that demonstrates a mech- 
anism by which magnetic fields tap rotational energy of a stellar core and expel the envelope. Our 
numerical setup, designed to focus on the basic physics of the outflow mechanism, consists of a solid, 
gravitating sphere, which may represent the compact core of a star, surrounded by an initially hydro- 
static envelope of ionized gas. The core is threaded by a dipolar magnetic field that also permeates 
the envelope. At the start of the simulation, the core begins to rotate at 10% of the escape speed. 
The magnetic field is sufficiently strong to drive a magneto-rotational explosion, whereby the entire 
envelope is expelled, confirming the expectation of analytical models. Furthermore, the dipolar nature 
of the field results in an explosion that is enhanced simultaneously along the rotation axis (a jet) and 
along the magnetic equator. While the initial condition is simplified, the simulation approximates 
circumstances that may arise in astrophysical objects such as Type II supernovae, gamma ray bursts, 
and proto-planetary nebulae. 

Subject headings: MHD — planetary nebulae: general — stars: evolution — stars: mass loss — stars: 
rotation — supernovae: general 



1. INTRODUCTION 

Over the last two decades magnetic fields have been 
identified as the principle, universal agent for creat- 
ing collimated astrophysical outflows. When a mag- 
netic field is anchored in a rapidly rotating object 
near the bottom of a gravitational potential, the field 
can act as a drive-belt, tapping rotational kinetic en- 
ergy and launching plasma back up the potential well. 
This magneto-rotational (MR) scenario for jet launch- 
ing has b een explored by numerous authors, both an- 
alytic ally (|Blandford fc Pavndll98l iPelletier fc Pudritzl 
Tflfll audTu~rnericallv JOuved Pudritz. fc Stonelll997t 



Krasnooolskv. Li. & Blandfordf 1999). These authors 
showed this mechanism can produce steady flows of mat- 
ter, in a sense that the outflow engine still operates while 
the fl ow is observable — as in j ets from young stellar ob- 
jects (iReipurth fc Ballvll200lf) and active galactic nuclei 
ijBegelman et al.lll984j) 

On the other hand, this mechanism can also operate in 
a transient event, linked to a rapid evolution of the source 
and driving an explosion. In this case, the engine loses a 
significant fraction of its power by the time the outflow 
(or its interaction with the environment) is detecte d. The 
propo sed mechanism for gamma ray bursts (GRBs: IPiranl 
2005) and supernovae (SNe), lie in the explosive regime . 
There is also growing evidence (B uiarrabal et al.ll200H) 
that collimated outflows from planetary nebulae (PNs) 
are explosive. 

Transient MR explosions have not been explored in 
as much detail as the steady-state models. In a col- 
lapsing star, differential rotation near the core may 
amplify t he field linearly (when turb ul ence is unim- 
porta nt; iKluzniak fc R.udermanl Il998t I Wheeler et"aT1 
2002) or exponentially (when turbulence is important; 

1 Dept. of Astronomy, U. of Virginia, Charlottesville VA, 22904; 
Levinson/VITA Fellow; seanmatt@virginia.edu 

2 Physics & Astronomy Dept., U. of Rochester, Rochester NY, 
14627; afrank@pas.rochester.edu, blackman@pas.rochester.edu 



Akivama etal] |2003; Blac kman et all |2006). When 
differential rotation twists a poloidal magnetic field, 
generating toroidal field B^, enough rotational en- 
ergy may be tapped to power a supernova explo- 
sion. The role of toroidal field pressure i n then driv- 
ing an outflow has been highlighted by iLvnden-Bell 
(1996), who explored magneto-static "magne tic tow- 
ers" to u nderstand jet propertie s (see also iLi et alJ 
I2001|) . and lUzdenskv fc McFadvenl l)2006j) has extended 
this work to explodin g stars . Numerical simulations 
bv iLeB lanc fc Wilson! (I1970T) and more recent ly by 
lArdelian. Bisnovatvi-Kogan. fc Moiseenkol l)2005f) sup- 
port the idea that MR explosions can be important for 
driving SNe and GRBs. 

However, the inherent time-dependence and general 
complexity of the system presents a significant challenge 
to our understanding of the basic mechanism and identifi- 
cation of key parameters. It is clear that a MR explosion 
is ultimately driven by the rotational kinetic energy ex- 
tracted from the material that is left behind (the stellar 
remnant). But, notably, there is still some uncertainty 
as to whether an accretion disk mus t form inside the star 
(as in lUzdenskv fc McFadver]l2006[) or wh ether the rota- 
tion o f the stellar core (e.g., a protopulsar: I Wheeler et alJ 
2002') alone will drive a MR explosion. 

In this Letter, we present a magnetohydrodynamic 
(MHD) simulation of a MR explosion. While our setup is 
simple, simil ar to the analyt i cal "p rotopulsar jet" model 
described bv I Wheeler et alJ l)2002jL our simulation cap- 
tures the nonlinear dynamics, as the magnetic field is 
twisted at the shear layer between the core and enve- 
lope. This heuristic approach enables a better under- 
standing of the basic MR physics and complements the 
more phenomenological approach of previous, more com- 
plex, numerical studies. 

2. PHYSICAL EXPERIMENT AND NUMERICAL METHOD 
2.1. Physics 
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The model consists of a gravitating, magnetized sphere 
(the "core," representing a stellar core) surrounded by an 
overlying envelope. The core is solid in a sense that no 
mass flows in or out of the surface, and it rotates as a 
solid body. The core is also a conductor with a rotation- 
axis-aligned dipole magnetic field anchored to the sur- 
face. This magnetic field couples to the envelope, whose 
behavior is characterized by the ideal MHD equations. 
This core-envelope model approximates the conditions 
in the interior of a star where there is a very steep den- 
sity gradient, for example in the proximity of a nearly 
degenerate core in an evolved star. There is an implicit 
assumption that the envelope self-gravity is negligible. 

We assume an initial state in which the envelope is not 
rotating and rests in hydrostatic equilibrium, and whose 
density falls of as r~ 2 . The initial dipole magnetic field is 
force-free and permeates the envelope. At a time defined 
as t = in the simulations, the core begins rotating at a 
constant rate of 10% of the escape speed. This initiates 
a shear at the boundary between the surface of the core 
and the base of the envelope. 

This strong initial shear and purely poloidal magnetic 
field approximates the conditions in a real star immedi- 
ately following an abrupt change in structure, such as 
a core collapse (or envelope expansion). If the collapse 
is rapid enough, subsequently amplified at the core- 
envelope interface will be much stronger than gen- 
erated elsewhere before the collapse, and the envelope 
rotation will be negligible. The envelope is unlikely to 
be at rest following such a change, and any bulk motion 
could quantitatively (but not qualitatively) change our 
results. For example, the ram pressure of an infalling en- 
velope would change the magnetic energy required for an 
explosion by approximately a factor of two. Finally, the 
field geometry is reasonable, since pulsars and magnetic 
white dwarfs are thought to have dipolar fields (though 
tilted with respect to their rotation axes). 

The poloidal magnetic field connecting the core to the 
envelope is twisted by the shear, generating an azimuthal 
component to the magnetic field, B^, near the inter- 
face, which increases approximately linearly in time. The 
magnetic pressure force associated with — Vi?^ will be 
directed generally outward from the core, so when B^ 
becomes dynamically important, the field will expand, 
driving the ionized envelope material in front of it. If 
the expansion is powerful enough, it can be explosive and 
drive off the entire envelope. We intentionally inhibit any 
sort of wind from the surface of the core, so the outflow 
in this explosion is composed of material that originally 
overlies the core (i.e., envelope material). The flow thus 
depends on the initial conditions existing above the sur- 
face of the core and is inherently a transient phenomenon. 
This contrasts with any steady-state wind theory whose 
solution depends only on the boundary conditions at the 
wind driving source. 

2.2. Parameters 

The key parameters in the system are the initial mag- 
netic, rotational kinetic, gravitational potential, and 
thermal energy densities near the core-envelope interface. 
It is instructive to compare the characteristic speeds as- 
sociated with these energies, namely the Alfven speed 
[va = B j '(47TP) 1 / 2 , where B is the dipole field strength 
at the equator of the core and p is the mass density at 



the base of the envelope] , rotation speed at the equator 
(u ro t), escape speed at the core surface (v csc ), and sound 
speed at the base of the envelope (c s ), respectively. 

The assumption of initial hydrostatic equilibrium cou- 
ples the thermal energy to the potential energy, effec- 
tively reducing the number of free parameters (by requir- 
ing the initial c s ~ 0.53w eS c)- Also, since the simulations 
can be scaled to any system with similar ratios of ener- 
gies, there are only two fundamental parameters, which 
can be cast as the dimensionless velocity ratios VA/v esc 
and v ro t/v csc . Here, we present results from a simulation 
with v A /v 0SC = 1.0 and v Iot /v esc = 0.1. 

The key parameter that determines whether or not the 
envelope will be driven off is x = (vAVmt) 1 ^ 2 /v csc - We 
have run a limited parameter study thus far (not pre- 
sented in this Letter), which suggests the threshold for 
explosion occurs near \ > 0.2. The model presented here 
has x ~ 0-3. 

2.3. Numerics 

We employ the 2.5-dimensional MHD code of iMattl 
(2002), an d the reader can fi nd d etails of the code ther e 
and also in lMatt et all <|2002fl and lMatt fc Balickl l|2004ft . 
The code, solves the ideal (non-resistive) MHD equations 
using a two-step Lax-Wend roff, finite difference scheme 
(Rich tmver Sz Morton! Il967) in cylindrical (137, (j>, z) ge- 
ometry. The formulation of the equations allows for a 
polytropic equation of state (we adopt 7 = 5/3), includes 
a source term in the momentum and energy equations for 
point-source gravity, and assumes axisymmetry (d/d<f> = 
for all quantities). 

The computational domain consists of five nested grids 
(or "boxes") in the cylindrical w-z plane. Each box con- 
tains 401 x400 gridpoints (in w and z, respectively) with 
constant grid spacing. The boxes are nested concentri- 
cally, so that the inner (first) box represents the smallest 
domain at the highest resolution. The next outer box 
represents twice the domain size with half the spatial res- 
olution (an so on for other boxes). The simulation ended 
before the explosion propagated beyond the fourth box. 
A circular boundary with a radius of 30.5 gridpoints and 
centered on the origin represents the surface of the core. 
Thus, the innermost box has a spatial resolution of 0.033 
core radii, R c , and a domain size of 13.1 i? c , and the 
fourth box has a resolution and domain size of 0.26 R c 
and 105 R c , respectively. 

We use standard outflow conditions on the outer box 
boundaries, appropriate for wind studies. We also use 
standard boundary conditions on the rotation axis (vd — 
0) and reflection symmetry on the equator (z = 0). The 
simulation domain, then, consists of a vj-z slice through 
a single quadrant. The surface of the core is repre- 
sented by a circular inner boundary, centered at the ori- 
gin (w = z = 0). Here, we enforce boundary conditions 
on a three gridcell layer, and we have performed several 
tests to ensure that this inner boundary behaves as ap- 
propriate for the surface of a solid, rotating, conducting 
and magnetized sphere. 

For the current study, we initialize the entire computa- 
tional domain with a spherically symmetric hydrostatic 
envelope, whose density falls off as r~ 2 (where r is the 
spherical radius). The envelope is initially stagnant (no 
rotation or motion) and is threaded everywhere by a 
dipole magnetic field that is anchored to the core, as 
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Fig. 1. — Greyscale images of log density (black is highest den- 
sity), poloidal field lines, and velocity vectors show the evolution of 
the system in the region near the core. The data are from the first 
simulation grid. The core is in the lower left of each panel, and the 
rotation axis is vertical. The four panels represent a time-sequence 
from 0, 1.5, 3.0, and 4.5 rotations of the core. 



Fig. 2. — Greyscale image of log density (black is highest density) 
far from the core, shown after 4.5 rotations of the core. The data 
are from the fourth simulation grid. The core is indicated by a 
white circle at the center. Vectors represent the flow velocity with 
the maximum vector length corresponding to 4.1 times v TO t- 



shown in the upper left panel of Figure ^ 

3. RESULTS: NATURE OF THE EXPLOSION 

At the start of the simulations, the dipolar magnetic 
field begins to wrap up. For the chosen parameters, 
becomes dynamically important after a few tenths of a 
rotation of the core. The panels of Figure ^ show the 
evolution during the first 4.5 rotations of the core. The 
expansion of newly generated is evident, as it drives 
the envelope material outward. This also expands the 
poloidal magnetic field (lines), though some field lines 
near the equator remain closed, where material is forced 
and held in corotation with the core. The region near the 
core is evacuated by the expansion of the magnetic field, 
indicative of a transient event, rather than the beginning 
of a steady wind. 

Figure [5] shows the density in the system after 4.5 ro- 
tations of the core, on a scale that is 10 times larger than 
in Figure n and the data have been reflected about the 
rotation axis and equator to better illustrate the flow. 
In the Figure, one sees a high-density swept-up shell of 
envelope material that is moving outward from the core 
and bounded by a shock. 

The shell has a quadrupolar shape. This symmetry is 
due to the coupling of the rotation to the dipole field, 
which produces a B^ that has a maximum strength at 
mid latitudes near the core, so magnetic pressure forces 
direct outward fr om the core, and als o toward the pole 
and equator (see lMatt fc Balicl3l2fl04D . This shape can 
be qualitatively understood as an explosion driven by 
detonating two doughnut-shaped charges placed at mid 
latitudes around the north and south hemispheres, which 
results in a convergent blast wave on the equator and 
pole. 



By 4.5 rotations of the core, the shell material near the 
rotation axis and equator is traveling faster than its own 
internal fast magnetosonic wave speed, as well as that 
of the ambient envelope ahead of the shell. We expect 
the expansion speed of the shell to roughly correlate with 
(v^va) 1 ^ 3 ■ At 4.5 rotations, the speed of the shell near 
the rotation axis roughly equals this. However, the flow 
speed of the shell material accelerates throughout the 
simulation, and does not seem to reach an asymptotic 
velocity by the end. The simulations were stopped after 
5.6 rotations of the core, since the steepening shock front 
began to produce numerical instabilities at that time. 

The total energy (gravitational potential plus magnetic 
plus kinetic plus thermal) integrated over the entire sim- 
ulation grid (excluding the core) increases monotonically 
in time and becomes positive at around 2.5 rotations. 
This means that there is enough energy for the entire 
envelope to escape from the gravitational potential well 
of the core. Indeed, by the end of the simulations, the 
shell material near the axis and equator is traveling faster 
than the local escape speed. On the other hand, mate- 
rial in the shell at mid latitudes is traveling substantially 
below its local escape speed. However, since this mid 
latitude material is still accelerating outward at the end 
of the simulation, and since the total energy in the sys- 
tem is well above zero, it appears likely that the entire 
envelope will escape from the core (if the simulation were 
able to run long enough). T his confirms the general pic- 
ture of lWheeler et alJ lf2002). In addition, the structure 
of the outflow suggest that, for lower x, a partial ejec- 
tion of the envelope may occur (a "failed" explosion) in 
which their is not only a jet (as in lWheeler et alJl2002j) 
but also an equatorial flow, leaving behind the material 
at mid latitudes. 

Finally, the magnetic field configuration after 4.5 rota- 



4 



Matt, Frank, & Blackman 




Fig. 3. — Three-dimensional rendering revealing the explosion 
mechanism, viewed from ~ 30 degrees above the magnetic equa- 
tor. The two blue surfaces are contours of constant density, each 
at the same density value. The dense shell of swept up envelope 
material (see Fig. |2J exists between the two surfaces. Gold wires 
trace the magnetic field lines and illustrate that the field is most 
highly twisted in the low density region interior to the shell. 



tions is apparent in Figure^ where we have exploited the 
axisymmetry of the system to create a 3-dimensional ren- 
dering of data from the fourth simulation grid. For scale, 
the outer isodensity surface corresponds to the outer edge 
of the shell in Figure El A comparison between Figures 
El and |31 reveals that the "hollow" region interior to the 
dense shell is actually filled with magnetic energy, mostly 
due to the component. The shell is like a magnetically 
inflated "balloon," which expands as the spinning core 
adds magnetic energy inside. Thus the spin energy of the 
core is ultimately being transferred to the kinetic energy 
of the envelope (and the core should spin down, as a re- 
sult). The interior region is magnetically (Poynting flux) 
dominated, as all material there flows more slowly than 
the fast magnetosonic speed. The flow along the axis is a 
magnetic to wer jet, similar to what is seen in othe r simu- 
lations fe.g.. lKato. Havashi. fc M atsumoto 2004), while 
the flow near the equatorial plane is more of a "magnetic 
pancake," a feature of the rotating dipole field and an 
interesting new result of this study. 

4. APPLICATION TO REAL STARS 

Our numerical setup is simple, which allows us to un- 
derstand the basic physical principles at work. The re- 
quirements for the mechanism to operate are a) a shear 
layer in which the inner region is rotating at a fraction 
of breakup speed and b) a magnetic field that results in 
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an Alfven speed comparable to the escape speed. Be- 
cause we focus on the MR mechanism, and ignore other 
physics (e.g., neutrino heating), we can scale our results 
to various astrophysical systems. This is instructive, as 
it determines what conditions would be necessary for a 
MR explosion alone to drive off the envelope. Here, we 
briefly apply our results to the progenitor stars of SNe 
and PNs. 

A simple model of a core-collapse SN suggests that core 
material conserves magnetic flux and angular momentum 
during collapse and naturally forms a rapidly-rotating, 
highly-magnetized proto-neutron star. We adopt a neu- 
tron star ("core") mass of 1.4 Mq and a core radius of 
R c = 10 6 cm. If the overlying envelope contains 1 Mq 
between r = R c and li?©, our simulation corresponds 
to a dipole field strength at the core surface equal to 
B = 3x 10 15 G. Note that the envelope mass is com- 
parable to the core mass, but the self-gravity of the en- 
velope near the core should not be important, and so 
the simulation still applies. In our simulation, we find 
that the core loses rotational energy at a constant rate 
of B 2 R^/3 = 3 x 10 48 ergs per rotation, which should 
extract most of the rotational energy of the core (~ 10 51 
ergs) in 1 second. These numbers are consistent with 
previous wor k and should be sufficient to eject the en- 
velope (e.g.. iWheeler et alJl2003 lArdelian et all 1200.1 
Blackman et aUKOOfH. 

The evolutionary phase preceding PN formation is 
marked by high mass loss rates, leading to an expan- 
sion of the stellar envelope, which overlies a proto-white 
dwarf. A shear layer at the cor e-envelope boundary is ex- 
pected ijBlackman et al.ll2001|) . and as the envelope den- 
sity decreases, the Alfven speed should rise. There may 
be a threshold density, below which the remaining en- 
velope may be expelled via a MR explosion. We adopt 
a white dwarf ("core") mass of 0.5 Mq, and R c = 10 9 
cm. At this time the envelope is essentially an extension 
of the massive outflow, and if it contains 0.05 Mq be- 
tween r = R c and 10 4 AU, our simulation corresponds to 
B — 9 x 10 6 G. Most of the rotational energy of the core 
(~ 10 48 ergs) should be transferred to the envelope in a 
few hundred years. The linear momentum of the swept 
up shell will be ~ 10 39 g cm s _1 , w hich is consistent 
with t he observations of proto-PNs by iBuiarrabal et alJ 
(2003). 

Finally, we note that the large-scale field in our sim- 
ulations is strong enough to suppress turbulence due to 
shear instabilities near the core, and the field therefore is 
neither subject to decay nor the exponential growth that 
is important in some other models. We have instead fo- 
cused on the physics of how the field actually drives the 
outflow. 
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